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The l/|k| contribution to the heavy quark potential is first generated at one 
loop order in QCD. We compute the two loop anomalous dimension for this 
| potential, and find that the renormalization group running is significant. The 

next-to-leading-log coefficient for the heavy quark production current near 
^r) ' threshold is determined. The velocity renormalization group result reproduces 

the a^ln 2 (a s ) "non-renormalization group logarithms" of Kniehl and Penin. 
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I. INTRODUCTION 



For systems involving a heavy quark-antiquark pair near threshold it is useful to combine 
the QCD coupling constant expansion with an expansion in powers of the relative velocity 
v. This expansion is facilitated by using non-relativistic QCD formulated as an effective 
field theory with an explicit power counting in v [PHI^H- For the potential between a heavy 
quark and antiquark this double expansion takes the form 

oo oo 

v{p,p r ) = v{n \ v(n) = E^ (n,i) > 

n=-2 j=l 

where V (n) ~ v n , V {n > ]) ~ v n a{ . (1) 

The even terms are first generated at tree-level (order a s ), and the odd terms \A 2fc+1 ) 
are first generated at one loop (order a 2 s ). 

Matrix elements for non-relativistic QCD systems typically depend on logarithms of the 
velocity v. For small v it is convenient to sum large logarithms of the form a s (mv)la(v) 
and a s (mv 2 ) ln(t>) by using renormalization group equations in the effective theory. This 
reorganizes the series in j in Eq. (||) so that: 

oo 

y{n) = y (n ' j) , where V (n ' j) ~ v 11 a s {m) j £ [a,(m) Hv)] k . (2) 

j k=0 

The simplest example of such a summation is the use of a running coupling constant, a s (/i), 
instead of the coupling at the matching scale a s (m). For fi < m the running coupling includes 
a series of a s (m) ln(/i/m) terms. However, it should be emphasized that the complete set of 
renormalization group logarithms are not determined by the simple replacement a s (m) — > 

A complication in summing the logarithms is the presence of two low energy scales: the 
typical momenta of the heavy quark ~ mv (soft scale), and typical energy ~ mv 2 (ultrasoft 
scale) . Two approaches have been proposed for dealing with the presence of two scales, one 
stage and two stage running. With two stage running, one matches QCD onto an effective 
theory at the scale fi = m and then runs to the scale fi = mv. At fi = mv one matches onto 
an effective theory called pNRQCD |8l,|T5| which has composite qq fields, and then considers 
the running to fi — mv 2 . In vNRQCD |14| , the running occurs in a single stage using a 



velocity renormalization group. The velocity renormalization group takes into account that 
the scales mv and mv 2 are tied together by the heavy quark equation of motion for all 
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H < m. We use dimensional regularization and the MS scheme with both the ultrasoft scale 
fiu ~ mv 2 and the soft scale fis ~ mv given in terms of a single subtraction point velocity v: 
fiu — 1711,2 and fis — mv - The renormalization group equations are written for the variable 
v. It is an interesting question as to whether both the one and two stage methods of running 
will sum the full set of a s ln(f) terms. In this paper only the single stage running will be 
considered. 

The running of the Coulomb potential at one and two loops, V^ -2,1 ) and \A~ 2 ' 2 ) ; is 
determined by the running coupling constant a s (fi) |]TB[ Q The summation of the leading 
logarithms for the the v potential was carried out in Refs. ]T8| , |iT| . In this paper we extend 
this analysis to the l/|k| potential by calculating its two loop anomalous dimension. 

If mv 2 ^> Aqcd, the v expansion can be applied to non-relativistic QCD systems in 
a perturbative manner. This is the case for ti production near threshold where mv 2 ~ 
4 GeV, and is the situation that will be analyzed in this paper. Of particular interest is the 
Coulombic regime where v ~ a s . In this regime the expansion in Eq. (JJ) has the form 



V 



y(-2,l)l + |y(-2,2) 



+ 



O'.s- 









+ 


V 2 


+ 



y(-2,3) + pK-1,2) + y(0,l) 
+ ... . 



+ 



4 + - + a.«° 
v v 



(3) 



If the a s ln(v) dependence is treated perturbatively these terms are referred to as the lead- 
ing order (LO), next-to-leading order (NLO), and next-to- next-to-leading order (NNLO) 
potentials. Note that since the l/|k| potential first occurs at one loop, it only contributes 
at NNLO. When the series in a s lnv are summed, the terms in Eq. (||]) will be referred to 
as leading- log (LL), next-to- leading log (NLL) and next-to-next-to-leading log (NNLL), re- 
spectively. In the Coulomb regime, the Coulomb potential must be kept to all orders. Each 
additional Coulomb insertion gives a a s /v 2 plus a factor of v (from the potential loop), so 
each new Coulomb interaction costs a factor of a s /v ~ 1. 

To study the threshold production of ti, a non-relativistic expansion must also be made 
for the electromagnetic production current^ 

C2 -Kp-p i xg-^WpV x g + .... 



2m 2 



(4) 



1 At three loops ultrasoft gluons can contribute to the running of the static potential |T7 
2 We will ignore effects associated with the top quark width. 
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The fields ifj* and x* create non-relativistic top quarks and antiquarks respectively. The c\ 
term contributes at order v°, and the C2 and C3 terms contribute at order v 2 . The current 
on the LHS of Eq. (H) is conserved, and has no anomalous dimension in QCD. However the 
non-relativistic current operators on the RHS are scale dependent in the effective theory, 
and the coefficients Cj therefore depend on logarithms of \i. The coefficients Cj each have an 
expansion in a s . The matching at fi ~ m is known to order a 2 for c\ pC| , pifl , and to order 
a s for C2 and C3 0, so the production current is known to NNLO with partial N 3 LO results. 
At LL order, one needs the tree-level matching C\ — 1 at fi — m, and the v 2 coefficients 
C2 and C3 can be set to zero. There is no one loop anomalous dimension for the operator 
ipt^X-pi so the LL result is that C\ — 1 at all v. 

At NLL order, we need the one loop matching for c\ at // = m, and the two loop running 
for C\. The coefficients 02,3 first enter at NNLL, at which order one would also need the 
three-loop anomalous dimension for C\. At two loops, the anomalous dimension for c\ was 
computed at the matching scale \i = m [f2C| . |2l|1 



9 ( \ 



-C F (±C F + ±C A y 2 s (m), (5) 



by studying the two loop matching condition for C\. For fi < m the anomalous dimension 
no longer has the simple form Eq. (|5]), but depends on the running of the quark potential. 
This anomalous dimension was computed in Ref. [O], and depends on the running values 
[y-(-2,i)]2 ) yi-^) x V^ - 1 ), and V^ 1 ' 2 ) (see Eq. below). It is interesting that the structure 
of this result implies that determining the RHS of Eq. (|) foifi<m requires the LL values 
of V ( - 2 ^ and V^, but the NLL value of V^ 2 \ The one loop running of V^ 2 ^ is well 
known, and the running of \/(~ 2 ' 2 ) was computed in Refs. |T8|JT9|1 . The two loop running 
of y* -1 ' 2 ) is computed in this paper. Using the running of these terms in the potential, 
we arrive at a complete NLL expression for c\. The running of the non-relativistic scalar 
current is also briefly discussed. 

In section II the effective theory is reviewed. We explain how reparameterization in- 
variance fixes the value of the lowest order coupling of an ultrasoft gluon to the Coulomb 
potential. The details of the computation of the two loop anomalous dimension for the l/|k| 
potential are given in section III, and in the appendices. In section III we give a derivation 
of the anomalous dimension using on-shell potentials, while in appendix [B| we repeat the 
derivation in the presence of off-shell potentials. Readers not interested in the technical 
details can skip to section IV, where our results are discussed. In section IV we expand 
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our renormalization group improved results in powers of a s to compare to finite order cal- 
culations in the literature. For the color singlet l/|k| potential the first a s ln(v) term in 
the series was computed in Ref. |15[, and our result for this term agrees with theirs. In 



Ref. ||22|| , Kniehl and Penin computed the a^ln 2 (a s ) terms in the wavefunction at the origin 
which they refer to as "non-renormalization group logarithms", since they do not involve 
factors of the /3-function for a s . We show that the second term in the series generated by 



our NLL production current agrees with the result in Ref. ||22|| . Thus, the solution of the 
renormalization group equations in the velocity renormalization group method does include 
these logarithms. Finally, we discuss our result for the NLL l/|k| potential and production 
current. 

II. THE vNRQCD LAGRANGIAN 

The vNRQCD effective Lagrangian has the form 

C = C u + C p + C s . (6) 

The ultrasoft Lagrangian C u involves the fields ip p which annihilate a quark, Xp which anni- 
hilate an antiquark, and which annihilate and create ultrasoft gluons. The potential La- 
grangian C p contains operators with four or more quark fields including the quark- antiquark 
potential. Finally, the soft Lagrangian C s contains all terms that involve soft particles which 
have energy and momenta of order mv. The terms we need in the ultrasoft Lagrangian in- 
clude 

The covariant derivative is = <9 M + igfifjA^ 1 = (D°, — D), so that D° = d° + ign^A , 
D = V — igfifjA, and involves only the ultrasoft gluon fields. The ultrasoft scale parameter 
Hu = mi' 2 , where v ~ v is the subtraction velocity. This v scaling for fijj is required for 
a consistent power counting in d dimensions |23]]. The covariant derivative on ip p and x P 
contain the color matrices T A and T A for the 3 and 3 representations, respectively. 

The Lagrangian £ p includes both the traditional quark potential and ultrasoft corrections 
to this potential which we will denote by C pu : 
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FIG. 1. QCD diagrams for tree level matching. 



£ P = - J2 V <*(^T (P, P) Ms ^p'l ^ P/3 X-P'A X-pr + £ 



(8) 



p,p 



The terms we need in C pu are fixed by reparameterization invariance |24| and will be de- 
scribed below. The on-shell potential V(p, p') a p\T has an expansion in a s and v, and 
a,/3,\,r denote color and spin indices. We will use the color basis in which the poten- 
tial V is written as a linear combination of 1 <8> 1 and T <g> T. The tree level diagrams in 
Fig. [I] generate terms of 0(v 2k a s ) in the QCD potential. The order v ~ 2 Coulomb potential 



is 



= (T A <g> T A )^- + (1 ® 1)^_ ; 
where the coefficients V^ 7 ' 1 ^ have an expansion in a s . The order v° potential includes 



(9) 



V<® = {T A ® T A ) 



vf } , vf ) (p 2 + p' 2 ) , vp a2 , y 



(T) 



+ 



+ (1® 1) 
where k = p' — p and 

_ (T 1 + (T 2 



v 2 w . v« 



+ 



A(p',p) 



2 m 2 k 2 
S 2 

.S-(p'xp) 



+ 



77T 



A(p',p) 



(T) 



77T 



T(k) 



(10) 



T(k) = o"! • o- 2 



3 k ■ (j\ k ■ er 2 



2 ,r ' k 2 x * k 2 

Matching the two diagrams in Fig. [I] to the v ~ 2 and v potentials at /x = m gives 

y( T ) = 4Wm) , Vi T) = Air a J m) , V (T) = 



ATTa s (m) 1 

1 7r a s (m) 

3 iV r 



V| T ^ = — 6ira s (m 
V« = , 



v 



(T) _ 7ra s (m) 



vf> = o. 



v« = o. 



(12) 



The LL values for the coefficients of the V^ 2 ^ and potentials will be needed below and 
are summarized in Appendix [A]. They are obtained by using the tree-level matching values 
in Eq. ([12]) and running using the one loop anomalous dimensions computed in Ref. Jl9| . 
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The order 1/v potential includes 

_2 



7T 



111 Ikl 



vf } (T A ®f A ) + V^(l® 1) 



(13) 



Tree level matching gives V^ l) = 0, and has zero one loop anomalous dimension |T3| 

so V^~^ = at LL order as well. At one loop in QCD, order 1/v potentials of the form in 
~~ |) are generated. The one loop matching of the on-shell potential at fi = m gives 



Eq 



23 



a 



2/ \ W 



(14) 



The color group theory factors are C F = (iV 2 - 1)/(2N C ), C A = N c , d = (iV c 2 - l)/(4iV 2 ), 
and Cd = N c — 4/N c . For the color singlet channel Eq. ( |T4"D agrees with Refs. p5| , |26|| . In 
the language of the threshold expansion JTOj the value of the coefficients in Eqs. (0) and 
(ITl) are from integrating out off-shell potential gluons at the hard scale \i ~ m where v — 1 
| |23| |. In the next section we will compute the two loop anomalous dimension for V^ -1 ), and 
determine the NLL value of the coefficients in Eq. (|i~3]). 

In addition we need ultrasoft corrections to the potential which are contained in C pu . 
Reparameterization invariance |24j] restricts the form of some of these terms by requiring 
that only the linear combination p — iD can appear. Here the covariant derivative acts on 
a quark or antiquark field with label p. The reparameterization invariant form of the T ®T 
Coulomb potential operator is 

rpA 



IP " P) 



1 

+ 2 



P ' - P + wy 

1>l T A iP P 



XL* T A X- P 

rpA 



(15) 



p' + p + iD) 2 

where D = ~3 + D . Terms in the v expansion are then generated by expanding Eq. (|15D 
with D <C p' - p. As written the ordering of color generators in Eq. ( |15"D is ambiguous. 
The correct ordering in the expansion is to write factors of D to the right of the T A , and 
factors of D to the left of the T A . Expanding Eq. (|l"5|) and keeping only the terms that we 
will need gives 



oj v( T ) f ABC 



pa 



k 4 



tk-V ^ 2 ^ 2l 



k 1 



v 2 +2 (k.v) 



2k 4 k 6 
zk • V V 2 
~k 4 2k 4 



+ 2 



k 6 



TV 



p • 



(16) 
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FIG. 2. Matching for the operator attaching an ultrasoft gluon to a potential interaction. 

where The first term couples an ultrasoft gluon to a four quark operator. 

The terms in the second and third lines in Eq. (JLQ) contain order 1/v and v° terms from 
the multipole expansion of the Coulomb potential. The first term could also have been 
determined from the on-shell matching calculation shown in Fig. |2|. However, by determining 
the terms in Eq. ( p!6|) using reparameterization invariance rather than matching we know 
that the coefficients remain equal to V c to all orders in perturbation theory. This saves us 
from the extra work that would be involved in computing the anomalous dimensions for 
these terms in C pu . 

The terms in the soft Lagrangian include 

£s = E { \fK - + <Pii<p* + e* <? c « } ( 17 ) 

-fas E {IVW-floJ^+IVW^ 1 ^ 
p,p' ,?,<?' 

+V [cv, c,]y (CT) ^ P + ( V T B V P ) (vVT"^)} + (V» - x, T - T ) . 

The fields and c g are the soft gluon and ghost fields, and tp q is a massless soft quark 
field with rif flavor components. U, W, Y, and Z are functions of (p, p', q, q') and matrices 
in spin and the index a denotes the relative order in the v expansion. For Feynman gauge 
and the case p 2 = p' 2 these functions were derived in Ref . [14 > |19[ . Beyond one loop terms 



proportional to (p 2 — p' 2 ) will be needed in C s , since besides its soft energy q° ~ mv the A q 
gluons can carry away a residual energy of order mv 2 . The LL values for the functions U, 
W, Y, and Z can be found in Appendix A of Ref. ||23|| . In addition, some NLL contributions 
to C s will be needed and will be discussed in section [111 A| . These additional contributions 
are obtained from one loop matching with two- loop renormalization group improvement. 

As an aside, note that it is not necessary to consider the ultrasoft renormalization of the 
vertices given in the soft Lagrangian in Eq. flT7p. One might think that diagrams such as 




(18) 



would effect the running of the coefficients in the soft Lagrangian. Diagrams analogous 
to the one in Eq. (|Hjj) , but with only soft gluons and quarks generate the running of the 
coefficient functions in the soft Lagrangian]] in Ref. | I9| . However, the graph in Eq. ([Tj|) has 



only one heavy quark, so a distinction between soft and ultrasoft gluons is unnecessary at 
this point. Noting that the soft vertices will always occur in pairs, it is in fact consistent to 
only dress pairs of the soft vertices by ultrasoft gluons: 






U " (19) 

Since these diagrams involve two heavy quarks, both types of gluons can occur. Since in the 
end it is only graphs such as Eq. ([19|) with two soft vertices that are relevant for constructing 
the theory, it is consistent to include the ultrasoft renormalization of soft vertices as a 
contribution to the four-quark operator in Eq. (|19D , rather than treating the subgraph as a 
contribution to the soft vertex, as in Eq. (|18D . Along with the diagrams in Eq. ([H]) there 
are graphs in which the ultrasoft gluon is exchanged between the two heavy quarks. 



3 In Ref. |L9[ the running of C s was calculated by examining loops with soft gluons that 
contribute to the Compton scattering process prior to integrating out the soft quarks. There 
it was noted that all ultraviolet divergent diagrams are in one-to-one correspondence with 
graphs in HQET, so that the running of this Lagrangian could be obtained from the known 
running in HQET |27 . 
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III. TWO LOOP ANOMALOUS DIMENSION FOR V^ 1] 



To calculate the NLL anomalous dimension for the l/|k| potentials we need to consider 
graphs in the effective theory of ordeif] aP s /v. These diagrams come in two classes, those 
with soft gluons, and those with a single ultrasoft gluon. The total anomalous dimension 
for the l/|k| potential is 7( 1,T ) = vdjdv V^ T \ Since fis = arid fiu = mu 2 , j^ 1 ^ can be 
written as the sum of a soft and ultrasoft anomalous dimension 

7^ = 7^ + 27^, (20) 

where Ts Hsd/dfis Vk ' an d 1u ^ud/d^u V*. ' • 111 ^ ne remainder of this section 
we discuss the computation of these two loop anomalous dimensions in detail. The calcu- 
lation is split into two parts, graphs with soft vertices and graphs with ultrasoft vertices. 
Graphs with soft gluons only contribute to 75, while those with an ultrasoft gluon contribute 
to both 7s and ju- 



A. Soft contributions 

The order ot\/v diagrams containing soft gluons that contribute to the anomalous di- 
mension are shown in Fig. |3]. The sum of diagrams forms a gauge invariant set. 

The two loop graphs in Figs. [3]a and [3]b involve an iteration of a potential and a soft loop. 
The vertices in these graphs are of LL order (tree level matching with one loop running) 
and are given in Ref. [J^]. In Figs. [5]a and [3|b, counting powers of v from the propagators 
and from the loop measures gives a if 1 , so the sum of powers of v for the three (amputated) 
vertices must give an overall 1/v 2 . The v scaling for the two soft vertices is o + a' — 2. 
In Fig. ^| we show two possibilities: a) has one V^~ 2 ^ insertion and two soft vertices from 
Eq. (0) such that a' + a = 2, and b) has one insertion of a V^°) potential from Eq. (]T0|) 
and two soft vertices such that a' = a = 0. We could also have a V^ 1 ^ potential plus two 
soft vertices where a + o' = 1; however this diagram is identically zero. The graphs where 
the potential and soft loop are exchanged simply give a factor of 2. 



4 This is the size (in v) of the amputated diagrams, so in contrast to the general power 
counting formula in Ref. JT1| we are not including the powers of v generated by external 
lines. 
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(a) 



(b) 










(c) (d) (e) 

FIG. 3. Order a^/v graphs containing soft gluons, ghosts, or massless quarks which are all 
denoted by a zigzag line. In (a) the dot denotes the Coulomb potential, while in (b) the dot 
denotes the order v° potential. The boxes denote soft vertices with insertions of the functions 
UW, yW, or Z^. In (a) the indices a + a' = 2. In (c), (d), and (e) the <g> vertex is 

obtained from the one loop matching in Fig. ||. 

The one loop graphs in Fig. |3|c, [5]d, and |3|e involve additional vertices in C s , denoted 
<8>, which are of NLL order (from one loop matching with two loop running). The one loop 
matching calculation for these vertices is sketched in Fig. f|. There are a large number of 
diagrams in the full theory (graphs on the left hand side), so only a few representative 
examples have been shown. To obtain the values for the operators on the right hand side 
we subtract purely soft effective theory diagrams from those in the full theory. To see how 
these operators arise, it is useful to recall that in the threshold expansion [fr|,|TT]| soft heavy 
quarks have a propagator 

P I _ i7r5 (q ) , (21) 



<?o + «e <?o 

where P stands for the principal value. In our approach, off-shell potential gluons and soft 
quarks are integrated out at the scale m when constructing the effective theory. When 
integrating out the soft heavy quarks the principal value term in Eq. (^1|) goes directly into 
a coefficient in the soft Lagrangian since this term is consistent with the scaling in the soft 
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J 3 3 



a) I +7 \ + 




b) A + + - 

^ % |_ 

FIG. 4. Contributions to the soft Lagrangian from one loop matching. In the full theory 
diagrams on the left the thick lines are massive quarks, while the thin lines are massless quarks. 
The matching in a) gives a s corrections to the functions and in Eq. (|i"7|). 

The matching in b) induces new operators that involve the scattering of soft gluons, ghosts, or 
light quarks off of a four fermion potential. 

regime, go ~ Tnv. For instance, in Eq. (|17D , Uqq = 1/go- The delta function contribution 
in Eq. ([H]) is associated with the potential regime since go ~ 0. When the delta function 
appears in a loop in the full theory (or threshold expansion) it forces gluons in the loop to 
have zero energy, or in other words to become potential gluons. It is these contributions 
which do not appear in the soft effective theory diagrams and must be made up by the 
operators shown on the right hand side of Fig. |j. 

The total contribution to the anomalous dimensions from the soft diagrams in Fig. [5] is 

m A)(7Ca — Cd) 3 , N 8CA(CA + Cd) 3/ , 
75 = 5 ot 8 (mv) a s [mv) , 

(i) A)Ci 3/ A WCaCx 3 
75= — 7T~ a s \ mv ) + a s \ mv ) , ( 22 ) 

where /3q = 11CU/3 — 4Tpnf/3, n/ is the number of massless soft quarks, and V c (^) = 
A7Ta s {mv) was used. In Eq. (p2|) the terms proportional to (3q can be inferred from Eq. (|14|). 
They simply turn the a s (m) , s in the matching result into running a s 's. At one and two 
loops terms proportional to the /3-function for a s completely determine the soft anomalous 
dimension for the Coulomb potential. However, the l/|k| potentials have additional contri- 
butions because the soft diagrams in Fig. [5]c, |3]d, and ^e have infrared divergences. The IR 
divergences from purely soft diagrams are not true IR divergences in the effective theory. 
For instance, in general they do not match up with IR divergences in QCD. The true IR 
divergences are from momenta < mv 2 , and the desired ultraviolet divergences are from mo- 
menta > m. Instead, the soft IR divergences are from momenta < mv, and match up with 



12 



ultrasoft UV divergences that are from momenta > mv to carry these UV divergences up to 

the hard scale. Writing the soft IR divergence 

1 1 / 1 1 \ 

— = , (23) 

cm e uv \£uv tins 

the first term contributes to the anomalous dimension in Eq. (p2j) . The 1/euv — V e iR term 

can be ignored since it simply takes the corresponding ultrasoft UV divergence up to the 

hard scale m|] Thus, to compute 75, all divergences from soft loops should be treated as UV 

divergences. The terms not proportional to (3q in Eq. (p2^ ) can be inferred from the result of 

the ultrasoft calculation in Eq. fl3~0|) of the next section. 

Despite the fact that Eq. ( p2]) can be inferred without a direct calculation it is worthwhile 

to examine the diagrams in Fig. |3|. Consider the graphs in Figs. |3]a and|3|b in Feynman gauge. 

The sub- loop with soft gluons is divergent, while the remaining potential loop is convergent. 

There is also a set of one loop diagrams (not shown) where the soft sub-loop is replaced 

by the one loop counterterms for V derived in Ref. |]19|| . We find that these counterterm 

graphs exactly cancel against a set of divergences in Fig. [|a and |3|b. For Fig. |3|b there is an 

exact cancellation, and so there is no operator mixing between and V^ l \ However, 

there are divergences that appear in Fig. [3]a that have no corresponding counterterm graphs. 

Consider the soft gluon case (the ghost and soft quark cases are similar). After performing 

the k° integration the loop integral for Fig. [3]a is 



1 



(E — q 2 /m)(p — q) 



d d t- 



jj(a) Jjfiu (a') 



(24) 



t 2 [(t°) 2 -(t + q-p') 2 ] 
where p and p' are the momenta of the incoming and outgoing quarks, E = p 2 /m and in 
Eq. ( f24"D the U's depend on q, t, p' and t°. For the divergences of interest performing the t 
integral (the soft loop) in d = 4 — 2e dimensions gives a factor of 

(p' 2 -q 2 ) 2 , s 

1 + . . . , (25) 



m 2 (p' - q) 4+2e e 



where the ellipsis denote order e° terms. The remaining loop integration is finite. For the 
one loop graph that corresponds to the soft sub-loop, the loop momenta q in Eq. (|25|) is 



5 This result, that there are no true IR divergences in the soft regime has not been proven 
to all orders in perturbation theory. However, its is likely that all IR divergences can be 



attributed to ultrasoft and collinear gluons in the spirit of the Coleman- Norton theorem p8 



plus IR divergences associated with the Coulomb regime that are reproduced by iterations 
of the potential. 



13 



replaced by p', and this diagram vanishes on-shell by energy conservation so no counterterm 
is generated. Performing the final integration and including the factor of 2 from the left-right 
symmetric graph gives 



Fig. 3a = i Vf) M c?M ^%(T A T B ®T A T B )\± + 2 In (A) + 



32m|k 



(26) 



The divergence in Eq. ( |26| ) contributes to the anomalous dimension for Vk- This seems 
slightly unusual because it was generated by a sub- divergence, whereas usually only the 
overall divergence in a diagram is relevant. In this diagram the loop integral with soft 
gluons generates a 1/e pole, and the remaining potential loop integral generates the l/|k| 
factor. The corresponding diagrams in QCD include graphs such as the vacuum polarization 
of one of the gluons in the box diagram. In this full theory graph the subdivergence due 
to the vacuum polarization insertion would be canceled by a counterterm. However, in the 
effective theory this divergence is instead absorbed into the coefficients of terms in the quark 
potential because the potential gluon components have been integrated out. Since matching 
the full theory box diagram gives a contribution to the l/|k| potential, it is not surprising 
that gluon vacuum polarization in the box graph contributes to the anomalous dimension 
of this potential. An alternative to the approach used to derive Eq. (p6|) is to use off-shell 
matching and running. In this case the soft anomalous dimension analysis will be different 
(since, for instance, the off-shell potential is gauge dependent). This approach is discussed in 
Appendix [B], where it is shown that the final answer for physical observables is unchanged. 

The diagrams in Figs. |3]c, |3|d, and |3|e are tedious to evaluate due to the large number 
of diagrams necessary for the matching calculation in Fig. [|.f] Since the graphs in Figs. |3]c, 
[|]d, and ||]e are one-loop diagrams there is no cancellation from counterterms. In Feynman 
gauge 

r 2,,2e „,3i 



Fig. 3c, d, e = i ^"^ (-CiT ' F n f {l <g> 1) + (C A T F n f - i C d T F n f ) (T A ® f A ) 
m\k\ Aire 13 v 3 ' 

+y C A C X {\ ® 1) - (I C A C d + g C\) (T A ®f A )Y (27) 



Only the terms proportional to nj in Eq. (|27|) have been checked by direct calculation. 
Also, in Eq. (p7|) we have assumed that the <S> operators simply run with the strong coupling 



6 In fact once a contribution to <S> has been identified it is simpler to directly evaluate the 
two-loop diagram obtained by combining the steps in Fig. ^ and Fig. |^. 
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FIG. 5. Order a^/v two loop graphs with an ultrasoft gluon and two potential insertions. 
The topologies shown give several different classes of diagrams depending on the vertices used as 
explained in the text. The diagrams obtained by flipping the graphs left-to-right and up-to-down 
are to be understood. 

constant, a s {mv). To motivate this, recall that a theory of propagating soft quarks and 
gluons has the same singularity structure as Heavy Quark Effective Theory (HQET). The 
operators that are generated in Fig. |4] correspond to diagrams in HQET with insertions of 
V 2 / (m) orap-A/m vertex. Since these operators do not run in HQET the operators that 
are constructed in Fig. |] should also simply run with the strong coupling constant. 



B. Ultrasoft contributions 



Next consider the order oe^/v diagrams with ultrasoft gluons. At one loop we have 
a single insertion of the V^~^ potential dressed by an ultrasoft gluon with gA° couplings. 
These diagrams are obviously zero in Coulomb gauge and in Feynman gauge it was shown in 



Ref. |14| that, with any potential, the sum of this set of one loop diagrams is also identically 
zero. At two loops in an arbitrary gauge there are many possible contributions. There 
are diagrams with two V^ 2 ^ vertices, and an ultrasoft tadpole generated by the seagull 
A 2 operator attached to one of the quark lines. These graphs do not have logarithmic 
divergences and therefore do not contribute to the anomalous dimension. There are also 
diagrams that are zero because they are odd in an ultrasoft momenta which we can omit. 
Next consider the topologies shown in Fig. [|. Several classes of diagrams are generated 
depending on the vertices used: 



1. vertex 1 and 2 are 2 \ and the ultrasoft gluon couples with p • A/ 



2. vertex 1 and 2 are 2 \ the ultrasoft gluon couples with gA°, and there are two 
insertions of p • V/m on quark lines, 

3. vertex 1 and 2 are V^~ 2 \ the ultrasoft gluon couples with gA°, and there is one insertion 
of V 2 /m on a quark line, 
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4. vertex 1 and 2 are V^~ 2 \ the ultrasoft gluon couples with g A , and there is one insertion 
of p 4 /m 3 on a quark line (together with the expansion of factors of the energy that 
appear in lower order diagrams which should be viewed as corrections to the effective 
theory states, see Refs. J7J^3|), 

5. vertex 1 is V^~ 2 \ vertex 2 is and the ultrasoft gluon couples with gA°, 

6. vertex 1 is the order 1/v potential from the multipole expansion of the Coulomb 
potential given in Eq. (|16D, vertex 2 is V^ 2 \ and the ultrasoft gluon couples with gA° 
vertices. 

7. vertex 1 is the order v° potential in Eq. (|i~6|), vertex 2 is V^~ 2 \ there is one insertion 
of p • V/m on a quark line, and the ultrasoft gluon couples with gA° vertices. 

8. vertex 1 and 2 are v° potentials from Eq. ([11]), and the ultrasoft gluon couples with 
gA° vertices. 

Insertions of operators with V's only need to be considered on quark propagators where the 
multipole expansion was used; together with the graphs in cases 6, 7, and 8 they build up 
the sub-leading terms in the multipole expansion. 

The graphs in cases 1, 2, and 5 do not give any contribution to the two loop anoma- 
lous dimension. The reason is that all ultraviolet divergences are exactly canceled by one 
loop counterterm graphs (these counterterms are generated by one loop graphs having one 
potential insertion dressed by an ultrasoft gluon, and in Feynman gauge were calculated 
in Refs. @,[19|). The graphs with the topology in Fig. ||]c,d are ultraviolet finite. For the 



remaining diagrams we can identify the corresponding counterterm graph by simply shrink- 
ing the smallest loop containing the ultrasoft gluon to a point. The only subtlety occurs in 
Fig. |5|h which has overlapping potential and ultrasoft loop integrals. For heavy scalars this 
graph was analyzed in Ref. |10| with the threshold expansion, while in an effective theory for 



heavy fermions this diagram was analyzed in Ref. |29[ for the case where the ultrasoft gluon 
is massive and couples with derivatives at the vertices. Since the nature of the overlapping 
integrals does not depend on the structure of the numerator or on having a massive gluon, 
this analysis will not be repeated here. In the effective theory this diagram comes with the 
appropriate factor of 2, so that it is canceled by the two one loop counterterm graphs that 
correspond to either making the loop integral on the right large or making the loop integral 
on the left large. Note that for a <ft 3 relativistic theory in d — 6 [30|, the sub-divergences 



for this diagram are also canceled in this way, but leave an overall divergence. For the 
non-relativistic effective theory diagram this overall divergence is not present. 
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For the graphs in case 4, the sum of insertions on quark lines inside the loop shared 
by the ultrasoft gluon are finite. The graphs with insertions on quark lines outside this 
loop are ultraviolet divergent, but are exactly canceled by the counterterm diagrams that 
correspond to shrinking the ultrasoft loop to a point. Therefore, the graphs in case 4 also 
do not contribute to the anomalous dimension. 

In a general gauge the graphs in cases 3, 6, 7, and 8 will contribute to the anomalous 
dimension. There are also additional graphs with the vertex in Eq. (|l|) which involves 
the coupling of an ultrasoft A* gluon to a potential vertex. Together this set of graphs 
gives a gauge invariant contribution to the anomalous dimension, so we can simplify their 
computation by choosing the most convenient gauge. We will choose Coulomb gauge since 
the graphs in cases 3, 6, 7, and 8 involve ultrasoft A gluons and vanish in this gauge. The 
result for the infinite part of the the remaining diagrams in Coulomb gauge is: 




, (28a) 



i aMfijyp? r i + A g fA mh) 

Sirm k e 4 




8tt m Ikl e 3 



(28c) 




8irm Ikl e 



2 -C A C x 1 ® 1 - ^f A {G A + C d )T A 



. (28d) 



Equations (|28|b,c) include the factor of four from their left-right and up-down mirror graphs. 
In the graphs in Eq. (EHf) the factor 



/4 e /4 e 



/4 e 



e|k| 1+2e ^ Ikl 



(29) 



occurs, where E = p 2 /m. The first term is included in the result displayed in Eq. (p8|), and 
the remaining terms are not shown. For v ~ v, we have [1$ ~ 1^1 ~ mv , A*c/ ~ E ~ mv 2 , 
so the soft and ultrasoft scale factors correctly minimize possible large logarithms in the 
effective theory. 

From the sum of diagrams in Eq. (^H|) we have the following result for the soft and 
ultrasoft anomalous dimensions of the order l/|k| potentials: 
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75 - 7c; - ^aIOa + C d J — , 

7s - lu ~ -^a^i • (3U) 

Recall that the soft and ultrasoft anomalous dimensions are given by differentiating with 
respect to In /is and ln/zjy respectively. 

It is interesting to ask how the result in Eq. (|30| ) would be reproduced if terms that vanish 
by the equations of motion were included in our effective Lagrangian. It is possible to include 
a term in our effective Lagrangian which makes the graph in Eq. (J28]b) give no contribution 
to the anomalous dimension, but the final result for observables remains invariant. This 
example is discussed in Appendix [B]. 

IV. RESULTS 

In this section the NLL heavy quark l/|k| potential and production current are discussed. 
After presenting the renormalization group improved results, we re-expand to compare to 
the finite order results in Refs. |15]| and ||22|| . We also discuss the behavior of the NLL l/|k| 



potential and the NLL production current as we run down from v = 1 to v = v, where v is 
the Coulombic velocity. 

A. The NLL l/|k| potentials 

The total soft and ultrasoft anomalous dimensions for the on-shell l/|k| potentials are 
obtained by adding Eqs. (^) and (p0|). Using the LL coefficient for the Coulomb potential, 
V c (v) = Aixctsimv), we find 

1 



lP = -—Po(7C A - C d ) [a s {mv)f - —C A (C A + C d ) [a s (mu)f 

4 



+ ^-C A (C A + C d ) a s (mv 2 )[a(mv)} 2 , 
An 

l ( s ] = Cx [a s {mv)f + - G A C X [a s {mv)f - - C A C X a s (mv 2 ) [a s (mu)} 2 , 

ATT 7T 7T 

lu ] = ^~Ca(C a + C d ) a s (mv 2 )[a s (mu)} 2 , 

T « = -- CUd a s (mu 2 )[a s (mu)] 2 . (31) 

7T 
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In full QCD (taking the QCD scale parameter /x = m), the first logarithm generated by 
the ultrasoft anomalous dimension corresponds to a \n(E/m), while the first logarithm 
generated by the soft anomalous dimension corresponds to a ln(|k|/m). For the velocity 
renormalization group, the total anomalous dimension is simply 75 + 27^: 



3tt 



C A (C A + C d ) [a s (mv)f 



(32) 



+ -C A (C A + C d ) a s {mv 2 )[a{mv)f 

7T 



1 



/3 Ci [a s (m^)f 



16 



C A Ci [a s (my)f 



24 



C A Ci a s (mz/ 2 )[a s (mi/)] 2 . 



Note that no other 1/v potentials are generated for v < 1 by operator mixing. Integrating 
Eq. ( |32"D using the one loop /3-function for a s and the one loop boundary condition in Eq. (|I~4"D 
gives 



(7C A - C d ) 



a 2 s (m) + 



{7C A -C d ) , 8C A (C A + C d ) 



3/3 



(^ 2 - 1) a 2 (m) 



(33) 



8Ca(Ca + C d 



where 



- 1 - 2 \n(w) a 2 s {m) 
( *i 16 C A d 



0o 



2 i\ 2( \ 4&C A Ci 2 



(^ - 1) a 2 (m) - 







2In(w) 



« s (mi/) 



a s (m^ 2 ) 



(34) 



a s (m) ' a s (mis) 2 — z 

Projecting onto the color singlet channel, = — CpV( T \ and setting C d = 8C F — 3C A 
and C\ = C A C F /2 — C 2 F gives 

8C A C F (C A + 2Cf) 



vl s \») = ( 



- C F C A ) a 2 (m) + 



8C A C F (C A + 2Cp) - 



0o 



r<2 

OfOj — — 

2 3/3 

1 — 2 ln(w) a 2 {m) . 



[z 2 - 1) a 2 (m) 



(35) 



The projection onto the color octet channel is vjf* = V^ X) + (C A /2 — Cp)^ 1 . 



Logarithmic corrections to the color singlet 1 - ) potential were also considered by Bram- 



billa et al. [13] using pNRQCD, but were not resummed. Brambilla et al. have 



-C F C A a 2 (r) 



4C A C F (C A + 2C F ) 
3tt 



a s (/x)a s (r) 2 ln(/xr) 



(36) 



To compare this expression to ours: a(r) — > a(mu), r — > l/{mv) and fi = fiu = m^ 2 since r 
corresponds to the soft scale and /x corresponds to the ultrasoft scale. Expanding Eq. 



m a s [m) gives 
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V { k s \v) = -C F C A a 2 s {m) + pC A a\{m) In[l/(mr)] 



7T 



AC aC F (C a + F ) 3 , . . 

a^(m) m(//r) + . . . . (37) 



3tt 

To compare this to our result we expand the resummed logarithms in Eq. (|35|) : 



VfV) = " C F C A y s {m) + a2(m) ln(„) - ^ a» ln(„) 

4CUC F (CU + 2C F ) _ , , 3 



« s (m) 3 ln(z/) + ... . (38) 

37T 

In Eq. (0), the first and second ln(z/) terms are entirely from the soft anomalous dimension 
in Eq. (|31~D , while the third ln(z/) term is from a combination of the ultrasoft and soft 
anomalous dimensions. The first ln(z/) term in Eq. ([38]) agrees with the ln[l/(rar)] = ln(z/) 
term in Eq. fl37|), and the third ln(^) term in Eq. ( p8[) agrees with the ln(//r) = ln(z/) term 
in Eq. (|37]). 

The second term in Eq. (|38| ) does not appear in Brambilla et al.'s expression in Eq. (j3"?D . 
This is because it depends on whether an on-shell or off-shell potential is used for the 
matching and running^]. We have used an on-shell potential, while Ref. [TJ|] uses off-shell 



Coulomb gauge which includes a = — 47ra s (r) potential as defined in Eq. (|B1| ) of 

Appendix B. The second logarithm in Eq. (^) should appear from the coefficient of this 
potential. Transforming Brambilla et al.'s potential to a potential using Eq. ( |B4| ) 
gives 

= ^fo? s {m) - al(m) ln[l/(ror)] + . . . . (39) 

2 Z7T 

The term with a ln[l/(mr)] = ln(z/) agrees with the second ln(z/) term in Eq. (|3"8"|). The sum 
of terms without logarithms in Eqs. Q57[ ) and (^) also agrees with Eq. fl3"8p. Note that the 
next term in the series in Eq. (|38D does not give agreement with Eqs. (|37D and (^) since 
Brambilla et al. did not attempt to systematically sum all the logarithms. 

To see the effect of the running on the value of the V^ -1 ) potential, consider the case 
of top quark production near threshold. Using a s {m t ) = 0.108 and Eq. (|35]), the one loop 
matching value is: 



7 The issue of on-shell versus off-shell potential was discussed in Ref. [23] , and in the context 
of the leading log summation is discussed further in Appendix 0. 
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FIG. 6. Values for the NLL running y( _1 ) potential for top quarks and n/ = 5 massless flavors. 
The solid and dashed lines are the color singlet and octet coefficients for the 7r 2 /(m|k|) potential. 
The solid vertical line marks the Coulombic regime where a s {mv) = v. 



v£°(i) 



-0.0362. 



(40) 



For a Coulombic system we determine the velocity v by solving a s (mv) = v. Using m t = 
175 GeV and the one loop running of a s ([i) with rif = 5 gives v = 0.145. The NLL color 
singlet and octet coefficients for the l/|k| potential are plotted in Fig. || At v — v the 
running coupling is 



-0.0027. 



(41) 



which is a substantial change from Eq. (|40f) . In comparison, the terms shown in Eq. 



where the combination a s ln(v) is treated perturbatively give V, 



-0.0313. Thus, 



the summation of logarithms for the l/|k| potential is quite important; it decreases the 
coefficient by an order of magnitude. Using the unexpanded results in Eqs. (|36|) and (39) 



gives V; 



(«)/ 



0.0201, so our resummed result is also quite different from the fixed order 



result in Ref. ill. 
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B. The NLL tt production current 



In terms of the running color singlet potentials the anomalous dimension for the produc- 



tion current in Eq. (|4]) is [p}] : 

d v^M (vj'Ku) , vW 



r y c Av) = v— ln[ci(z/)l - , , 

,clK J dv L V ;J 16tt 2 \ 4 



(42) 



For the vector production current in Eq. (|]) only spin-1 states are produced, and we can set 
S 2 = S(S + 1) = 2. However, our analysis also applies to the scalar production current 



c 0) Z ^x*- 



p • 



(43) 



which, for example, contributes to the process 77 — > tt [|3T]. Therefore, we will keep the 
factors of S 2 explicit in our results. The boundary condition is given by the matching 
condition at v = 1 (/i = m) which is known to two loops plfl . For the NLL approximation 
only the one loop matching condition should be used. For the vector current 

2CVa s (m) 



d(l) = 1 - 



7T 



(44) 



Integrating Eq. ( [42]) with the running potentials in Appendix [A], Eq. ( |A1| ) and in Eq. (|33] 
we find 



In 



d(l) 



di 7ra s (m)y lj + a 2 ira s (m)(l — z) + a 3 ira s (m) \n(z) 



+0,4 ita s (m) 
+a 6 ixa s (m) 
+a 7 Tia s {m) 



1 _ ^1-13^/(6/30) 
2 1 



+ 05 7ra s (m) 



1 _ z l-2C A /Po 
1 



^1 _ i i n 2 (2) _ l n ( w ) In f-^-) - Li 2 f 

W ln(w) - -ln(2w- L 
2w — 1 2 



2w 



(45) 



where z = a s (mv) / a s (m) and w = a s (mv 2 ) ja s {mv\ The coefficients aj in Eq. ( f45| ) are 
32CyiCp'(C J 4 + ICp) 



a 1 



«2 



a 3 



3/3 2 

-C F [3/? (26Cl + 19C A Cjr - 32C 2 F ) + C A (208^1 + 597C A C F + 716Cf.)] 



3/3 2 (6/5 -13C A )(/?o-2C A ) 



78/? 2 C A 

2C2(66/5o-179C A )(/3 -2C A ) 
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-CaC f 



6(49 - 3S 2 )/3 2 - (1126 - 39S 2 )(3 C A + 1067C A 



-24C 2 A (6f3 -UC A )(f3 -2C / 



-24C 2 (3f3 -llC A )(5C A + 8C F ) 



13C^(6/3 -13C A ) 2 



a 5 



C 2 



(4S 2 - 3)A) + (15 - 14S 2 )C A 



6(A) - 2C A ) 2 



a 6 



-16C 2 (C A + 2C F ) 



A further check on our result can be made by comparing it to the a^\n 2 (a s ) corrections 
calculated by Kniehl and Penin in Ref. [22|. Near threshold the ti cross section depends on 
the product |3"2]] : 



16C A CV(C A + 2C F ) 



(46) 



acx\ Cl \ 2 G c (0,0,E) \^(0)\ 



(47) 



where ip c is the leading order Coulomb wavefunction, and Gc is the Coulomb Green's 
function. In the approach used in Ref. p2[ , Gc embodies corrections to the wavefunction at 
the origin, (^(O)! 2 = I^^CO) | 2 [1 + A-^ 2 (0)] , and includes the large logarithms. They calculate 
the a 2 ln(a s ) and aj?ln 2 (a s ) terms and find: 

2 



A^ 2 (0) = -C F aj\n(a s 



2--S 2 

3 



Cp + C j. 



c 



^aUn 2 (a s )\^C 2 + 



71 



41 



7 



— - — S 2 
12 12 b 



CfCa + 



(48) 



for the terms not involving /3q. In our approach the large logarithms in the cross section 
all appear in the running coefficient C\(v), so we expect that the logarithms of Kniehl and 
Penin will be reproduced by 

2 



A^ 2 (0) 



cx 1 



2 In(a.) 7ci (l) + ln 2 (a s ){ 7 ;(l) + 2[ 7ci (l)] 2 } + . . . 



(49) 



where we have expanded to second order in ln(z/) = ln(a s ). The [7ci(l)] 2 term does not 
contribute at order c^ln 2 a s , and can be dropped. The remaining terms in Eq. ([49]) give: 

2. 



A^ 2 (0) = -C F a 2 s \n(a 



2--S 2 

3 



(50) 



C 

TV 



^ln 2 (a s )<|^ 2 + 



3 ( 



41 7 

|_g2 

12 12 



2 
-( 

3 



CfCa + -C\ — ^ 



2 



2--S 2 )( < 



which agrees exactly with the result from Ref. [^] in Eq. (48") after setting f3 = 0. Thus, 
we have shown that the logarithms of Kniehl and Penin can indeed be associated with 
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FIG. 7. Values for the tt production current coupling c\{u) for nj = 5 massless flavors. At 
LO and LL c\ = 1. The large and small dashed lines show the value of c\ at NLO and NNLO 



[21] respectively. The solid line shows the running value at NLL order given in Eq. (45). The 
dotted line shows the result of promoting the coupling in the NLO result to a running coupling 
(the approximation used in Ref. |^3|). The solid vertical line marks the Coulombic regime where 
a s (mu) = v. 

renormalization group logarithms. Expanding Eq. (f4"5|) to higher orders gives the a^ln 3 (a s ), 
al ln 4 (a s ), etc. terms. 

For QCD with n/ = 5, the coefficients in Eq. (pi) are 

4.308 - 0.417 S 2 , 



«5 



4.113 



«2 



-2.173, 



a 3 



a 4 = 5.731 



2.347- 1.209 



a 6 = -0.914. 



«7 



6.170 



(51) 



The running coefficient C\(v) is shown in Fig. [7] for S 2 = 2. For comparison we have 
also shown by a dotted line the value of the running coefficient obtained in Ref. [|33] by 
an approach based on simply taking a s — > a s (mv) to approximate the NLL value for the 
coupling. This approximation misses many of the a s ln(t> ) terms, and does not provide a 
good estimate of the NLL result, as seen in Fig. [7[ Some of the more important logarithms in 
our NLL result include the large Xwimv 2 jm) terms that enter through the mixing generated 
by ultrasoft gluon diagrams in the potential. 

We have also shown the NNLO value for c\ in Fig. [7[ As pointed out in Ref. |2l| the 
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NNLO value of ci(l) is fairly large. From Fig. |7|we see that summing the logarithms reduced 
the size of the NLO matching correction by a factor of two. It would be interesting to see 
if the running induced at NNLL continues to improve the convergence of the expansion. A 
consistent calculation at this order requires the three- loop anomalous dimension of C\, as 
well as one loop running of the v 2 coefficients C2 and C3. This computation appears quite 
involved since running of the current at this order will likely depend on the running of higher 
order terms in the potential. 

Several groups have analyzed the ti cross section predictions at NNLO using effective field 
theory techniques [p3 H37|j3T|| . It should be straightforward to incorporate the renormalization 
group improved current and potentials into their analysis by simply choosing a value of v 
appropriate to the threshold region. The value of v only needs to be of order the velocities 
in this region for the large logarithms to be minimized. Additional logarithms that appear 
in evaluating matrix elements with the potential will not involve large ratios of scales since 
they are of the form ln(E/nu) and ln(|k|/jUg) where \x\j = mv 2 and [i$ = mv. In a more 
precise analysis one might wish to use a s (mv)C F = v to determine the value v to use. 

We would like to thank I. Rothstein and J. Soto for discussions. This work was sup- 
ported in part by the Department of Energy under grant DOE-FG03-97ER40546, and by 
the National Science Foundation under NYI award PHY-9457911. 



APPENDIX A: LL values for the order v 2 and v° potentials 



The leading log values for the coefficients of the 2 -* and potentials are [19|: 



yf)(z,) =A7ra s {m)z, 
V^ T \u) = 47r a s (m)z —^-a s {m) 



(Al) 



3/3o 



l-z 



64nC A 
Wo 



a s (m) In (to) 



71 



C A {352C F + 91C d - 144CU) - 3/3 (33C A + 32C F ) 



+ 



+ 



39A,C A 

8vr(3/? - HC A )(5C A + 8C F )a s (m) 
ISC A (6(3 - ISC A ) 



a s {m) 



z-l 



7t(/3 - 5C A )a s (m) 



(Po ~ 2C A 

28TTd 



Z (\-2C A /P ) _ 1 



^(1-13^/(6/30)) _ 1 

8-K(AC F +C d -3C A ) 



a Am 







) In {w^j 



3/3o 



a s (m)( 1 — z) + 



a s (m) In (w^j 
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vfV) 
vi 1} (^) 



27ra s (m) 
(2C A - ft) 

7ra,(m), 

c 

_vr« £ (m) (1 _ 2Ca/A)) 
3 

27ra s ('m > ) 



+ — vra s (m) 



z _ 4 z (i-c A /Po) 



where z = a s (mv)/a(m), w = a s {mv 2 ) / a{mv), and we have included the N c dependent 
terms that come from matching the tree level annihilation diagrams. In the color singlet 
channel, Vfiy) and Va(^) were first calculated in Ref. [18] and agree with Eq. (|A1|). 



APPENDIX B: POTENTIAL OPERATORS THAT VANISH ON-SHELL 

It is interesting to consider how the soft and ultrasoft anomalous dimensions for V^ 1 ^ 
are affected when operators that vanish on-shell are included in the potential Lagrangian, 
C p . These new operators can have non-trivial anomalous dimensions. In this Appendix we 
consider two examples of how including these operators affects intermediate results, but in 
the end yield the same result as the on-shell potential for observables. Equivalently, it is 
shown that running the Lagrangian with the new operators from v — 1 to v — v and then 
removing them by field redefinitions or operator identities reproduces the on-shell running 
value of V^l 

Consider the scattering Q(pi,p) + Qip®, ~ p) — > QG°3>p') + Q(pI,~ p')- As our first 
example, consider including in the v potential in Eq. ( |10D terms of the form: 



y(o) 



+ 



(pg-pg) 2 

k 4 



V A V (T ® T ) + Vti (1 



4m 2 k 4 



:bi) 



We will refer to this potential as an off-shell potential, since on-shell it vanishes by energy 
conservation, where p® = pi, p' 2 = p 2 . Define Va = Vj^' 1 ^ + vjQ . In Ref. |23j it was 



shown that there is an operator identity whereby the time ordered product of a V c and a 
Va potential gives Vk potentials: 



(B2) 



\/ = " c — T A T B ®T A T B + c T A ® T A + . . . . 

32mA; 32mA; 
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where the ellipses denote additional terms which vanish on-shell. The identity in Eq. ( |B2|) 
allows us to remove the potential in Eq. (|ET| ) at an arbitrary velocity scale v in favor of a 
purely on-shell potential. Transforming iy) — ¥ gives 



V, 



en, 



v, 



V) 



vfV) + 



32tt 2 
1 

32^2 



(B3) 

induces a color singlet l/|k| 



In a situation where (v) 
potential of the form 



transforming 



C 2 



32tt 2 c 



u)-^V^\v)V { P{p). (B4) 
'' potentials were also considered in 



v 



(Transformations between order 1/v and 
Refs. |3|J3|J3|].) 

To make the implications of Eq. ( B3|) more clear, we will consider an example and show 

(T 1) 

that including the ' potential will not effect predictions for observables when running 
below m. Consider matching offshell in Feynman gauge where the matching coefficients at 
v — 1 are = 4:iia s (m) and =0, and the one loop anomalous dimensions for the 
potentials in Eq. (Bl|) are p3 



-2p a s (mu) 







(i) 



0. 



(B5) 



This anomalous dimension has contributions from the one loop graphs with two soft gluons, 
ghosts or quarks where: the two factors of (p /2 — p 2 ) are from the soft vertices, or these 
factors are from two insertions on the soft propagators, or where one factor is from the 
propagator and one from the vertex (for the ghost loop). The solution of Eq. (|B5|) is 



(T), 



V) 



and Vl 1} (z/) = 0. The one loop counterterm which generates this 



running affects our soft anomalous dimension computation since now there are counterterms 
of the form 



(P 



/ 2 



P 



2\2 



m 2 (p' 



P ) 4 6 



(B6) 



These counterterm give rise to one loop diagrams which exactly cancel the divergences in 
Eq. ( ^6|) for the two loop graph in Fig. §a. Recalling that the one loop matching values for 
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the l/|k| potentials are also different [23|, we find that the running V^ T (v) coefficients in 
Eq. (p3[) become 



V, 



(T), 



V) 



(3C A - C d ) 



a 2 s {m) + 



(3C A -C d ) , 8C A {C A + C d ) 



8C A (C A + C d 



0o 



4 3[3 
2 \a(w) a 2 (m) , 



(z 2 - 1) a 2 (m) 



(B7) 



Vl 1 \v) = C 1 at l a (m) + 



Ci- 



16 C A d 



A) 



[z 2 - 1) a 2 (m) 



48C A Ci 



a 2 (m) z — 1 — 2 ln(w) 



This Feynman gauge off-shell potential is consistent with the transformation in Eq. ( p3[ 

47ra s (mi/) and Eq. ([B7|) , into the on-shell case 



(T) 

which transforms the off-shell case, V" A 



V, 



and Eq. (0). The comparison in section |IV A| with Ref. [JT5j| shows that our 



on-shell analysis is also consistent with an off-shell Coulomb gauge potential. 

With the off-shell potential in Eq. ( |B1| ) there are two new graphs which contribute to 
the anomalous dimension for the production current: 

fB8) 





and the anomalous dimension is therefore 



d 



-V«, 



V) 



16vr 2 



64vr 2 



(B9) 



where V 



(*)/ 



V^V) - C F V. 



Since the last two terms in Eq. ( P^| ) are invariant 
under the transformation in Eq. ( |B4| ) the off-shell potential gives the same prediction for 
the running of the production current. In calculating observables the operator identity in 
Eq. (|B2| ) guarantees that the time ordered product of the running Va an d produces 
the same effect as the soft contribution to the running of the on-shell in Eq. fl26|) . 

As our second example, consider including an operator which would vanish by the lowest 
order free equations of motion: 

Vf 



E 



p,p' 



p' - p) 



4, T A (id ~^)% xl p/ f A X - P + X <- V ■ 



(BIO) 



The only feature of the operator in Eq. ( |B10| ) that is different from typical operators that 
vanish by the equations of motion (for example, in HQET P7|| ) is its non-local nature relative 
to the scales p ~ p' ~ mv. At one loop there are ultrasoft gluon graphs which mix into the 



operator in Eq. ( [B 1 0|) , for instance from the diagrams: 
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In the first graph the vertex is the first term in Eq. ([Tq) and the ultrasoft gluon is A, while 
in the second diagram the cross is a insertion of the ultrasoft V 2 /m operator in Eq. (|7|) and 
the gluon is an A . In Feynman gauge both graphs produce a divergence of the form 

l A=^4. (B») 
e (p' - p) 2 

which vanishes by the equations of motion, but off-shell renormalizes the coefficient of the 
operator in Eq. flB10|) . To see how this effects the calculation of our ultrasoft anomalous 
dimension, consider Coulomb gauge, where only Eq. ( |Bll| a) is non-zero. In this case there 
is a one loop counterterm diagram which exactly cancels the divergence in Fig. (|28|b). After 
running down to the low scale the value of ' (v) are different, but our potential also 
includes the operator in Eq. (P10|) with coefficient Vf{v)- If we wish to remove the operator 
in Eq. ( |B10| ) we can do so with a field redefinition 

i>l - ^ + E l^^A-T^-v" fA X~ P , (B13) 

which induces a term from the ipUido — p 2 /(2m)]^ p Lagrangian which cancels Eq. (|B10[) . 
This field redefinition also gives other new contributions. For us the important point is that 
the Coulomb potential induces a six-quark operator of the form: 

" £ 5 ♦^V^<«^V_^_ p)! • (B14) 

(The contraction of color indices and factors of T A have been suppressed.) Usually a six 
quark operator could not possibly effect the running of a four quark operator such as the 1/ |k| 
potential. However, because of the momentum dependence in the denominator of Eq. ( |B14|) 
it induces a non-zero tadpole diagram where the fields X-p' an d X-p> are contracted. This 
tadpole graph produces a l/|k| and makes up for the running in 14 that was removed when 
the contribution from Eq. (|28|b) was canceled by a VV counterterm. 
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